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are introduced by making one tunneling rate energy-dependent. We show that the non-Markovian 
master equation may fail if these effects are strong. We evaluate the stationary current, the zero 
frequency current noise and the occupation dynamics of the resonant level by means of a quantum 
master equation approach within different approximation schemes and compare the results to the 
exact solution obtained by scattering theory and Green's functions. 
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I. INTRODUCTION 

The ongoing progress in measuring tiny fluctuations of 
charge currents through nanoscale conductors^! 3 - has led 
to an increased theoretical interest in non-Markovian ef- 
fects revealed in such experiments 4 -^. Of particular inter- 
est are the non-Markovian dynamics induced by the cou- 
pling to fermionio 6 -! 7 - and/or bosonio^i 1 ^ environments 
and their influence on steady-state transport observables 
such as the cumulants of the stochastic charge transfer 
process^. Our present work focuses on such effects due 
to the coupling to electronic reservoirs. 

There exist various techniques for describing open 
quantum systems coupled to fermionic reservoirs: e.g., 
scattering theory and Green's function o 12 ' 13 ' 14 , or quan- 
tum master equations starting from the von-Neumann 
equation for the total density operator—, or the 
Wigncr-Boltzmann approac h 16 ! 17 ' 18 . Master equations 
are widely considered with a conductor-lead coupling 
in Born-Markov approximation (see e.g i 19 ' 20 ' 21 ), which is 
only strictly valid for weak coupling and constant contact 
density of states in the energy range of interest. Con- 
sequently, to explore effects beyond the common Born- 
Markov-approximation, one would study higher-order 
perturbation theory in the contact coupling (e.g . 22 ' 23 ' 24 ) 
and/or allow for energy-dependent tunneling rates, i.e., 
go beyond the wide-band approximatio n 14 ' 25 ' 26 ' 27 ' 28 ' 29 . 
In this work we choose the latter whilst analyzing pertur- 
bative approaches in lowest-order tunnel coupling only. 
In particular, we address the question of in which lim- 
its it is possible to describe non-Markovian physics with 
master equations. 

A conceptionally simple model for this purpose con- 
sists of a single resonant level (e.g. the ground state 
of a quantum dot or a molecule) weakly coupled to two 
electronic leads in equilibrium and a Lorentzian-shaped 
density of states for one of the leads. This model is equiv- 
alent to two serially coupled quantum dots^ and can be 
interpreted as a quantum dot coupled to a reservoir with 
finite electron relaxation time. 

To enable an exact solution, and for the sake of simplic- 
ity, we neglect Coulomb interaction and, consequently, 
prominent effects like Coulomb blockade or Kondo cor- 
relations. This, of course, constrains the use of our re- 



sults for quantitative understandings of transport exper- 
iments. Our aim, however, is a comparison of different 
approximation schemes for master equations in an elec- 
tronic transport problem that has an exact solution. 

We will show that the non-Markovian master equa- 
tion (NMME) in the wide-band limit produces reason- 
able results for the current and noise. However, reducing 
the band-width yields qualitative and quantitative devi- 
ations — even the emergence of unphysical results, like 
negative Fano factors. The exact time evolution of the 
resonant level occupation can only be obtained by the 
NMME in the wide-band limit where even a Markovian 
master equation covers the exact dynamics. In the short- 
time limit, the NMME result well approximates the exact 
evolution regardless of the bandwith. However, for very 
small bandwidths, the NMME generates negative, un- 
physical probabilities. We demonstrate how this can be 
avoided by a dynamical coarse graining method. 

The paper is organized as follows: In Sees. II. A. and 
II. B. we introduce the model and provide the known ex- 
act solution obtained by scattering theory and Green's 
functions. In Sec. II. C. the equivalence to the double- 
dot model is discussed. In Sec. III. A. we introduce 
the non-Markovian master equation and in Sec. III.B. 
the dynamical-coarse-graining approach. In Sec. IV the 
steady-state current and the Fano factor are compared to 
the exact and NMME solution. Finally, the occupation 
dynamics of the resonant level is discussed. 

II. MODEL 

A. Hamiltonian 

We start from the well-known single resonant level 
model^i 3 ^ 3 ^ 3 ^ 3 -! 3 ^ that is described by a Hamiltonian 
where tunneling between two leads can be realized via a 
localized quantum dot state, 

H = e d d ] d + ^2 tkac\ a c ka + ^2 (tkacftcka + h.c.) . (1) 

k,a k,a 

The fermionic operators d/dJ annihilate/create an elec- 
tron on the dot and Cfc a /ct annihilate/create an elec- 
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FIG. 1: Single resonant level coupled to two leads with one 
constant and one Lorentzian tunneling rate. In the infinite 
bias limit, hl — ► oo and /j,r — > — oo. 



tron in lead a £ {R,L} with momentum k. We choose a 
Lorentzian-shaped right tunneling rate, like i n 25 i 26 : 



r 2 



£fl) 2 + <S 2 



(2) 



and a flat left tunneling rate T^, both re- 
lated to the microscopic parameters via 
T a (uj) = 2irJ2k\ t ka\ 2 8(u - e ka )- A sketch of the 
model is given in Fig. Q] We work in the infinite 
bias limit, where we couple the left and right lead to 
two particle reservoirs in a way that the left lead is 
completely occupied and the right lead is completely 
empty, i.e., //l - ► oo and fj,R — > — oo such that the left 
Fermi function takes the value 1 and the right one takes 
the value 0. 

The retarded and advanced self energies of our model 
are evaluated exactl y 14 ' 32 i 34 as 



\tfca 



k,a 



UJ 



e ka ± i0+ 



A(w) =f i 



T(uj) 



(3) 



The imaginary part T(u>) = Tl(uj) + Tr(u>) is connected 
to the inverse lifetime, whereas the real part 



AH 



1 

2^' 



du) 



(4) 



(where V denotes the principal value) induces a level shift 
which is fully determined by the tunneling rate (Kramers- 
Kronig relation). It will turn out that all quantities con- 
cerning transport statistics and the occupation of the 
quantum dot can be expressed in terms of T a (uj) and 
e d only. 



B. Exact solution 

To evaluate transport quantities we need the transmis- 
sion coefficient 3 -^ , 



T( W ) 



|u;-e d -E*(a,)| 2 ' 



(5) 



which is related to the dot spectral function 
A(uj) := i [G r (uj) - G a (oj)] via T(u) = ^^f^A(uj). 

Scattering theory^ then yields the current I and the 
zero- frequency noise S which at infinite bias are 



oo 



s 



2^ 



TH[i-TH], 



(6) 



where e is the elementary charge (chosen negative). The 
time-dependent occupation probability rid(t) can be ex- 
pressed with Green's function a 36 i 37 . We use that at infi- 
nite bias the left Fermi function is unity and the right is 
zero. For simplicity we assume n<j(0) = and thus 



k 



(7) 



where we have used G^ ad {-t) = G R M (t) . 

The required Green's functions are obtained using 
equations of motio n 14 i 32 i 33 ' 34 and read 



G d,kL(t) - tkL 



du> 



1 



2-7T to — EkL + i0 + oj — e d — S fl (w) ' 



(8) 

which, upon inserting 1 = J_ 8(u> — e^i^dui, yields the 
explicit result 



n d (t) 



du , 



Tr X 



oc 

/ 



2tt 
du' 



-ioj t 



2ir uj' - uj + i0+ w'-e d - £ R (w') 



(9) 



This can alternatively be obtained from a direct calcula- 
tion without Green's functions 3 -^ 

For the Lorentzian-shaped right tunneling rate (|2|) , the 
Kramers-Kronig-relation (U|) yields the level-shift func- 
tion 



AM 



r 2 

1 Rfi 



UJ - S R 



2tt {uj-e R y+5l' 
and the self energy ^ can be simplified to 

r L , r 2 i 



£*M(o,) = Tl 



2-k UJ — £R ± iSr 



(10) 



(11) 



With Lorentzian shaped tunneling rates we can analyt- 
ically integrate the expressions for current and noise, 



3 



((6]), and with the abbreviations r := Tl + 26 R and 

e := ed — e R , wc obtain 



I = e 



2T L S R T T 2 R0 /ir 



(T L S R + r 2 R0 /7T)T 2 + Ae 2 T L S R ' 



g = /efl- 2W ^x 
I Trr 



4e 2 (r| + 8^) + f 3 (T 2 + 2r L J fl + 2T\^k) 

[(r L 6 R + r 2 R0 /n)r 2 + 4e 2 r L s R ] 2 



(12) 



C. Comparison with double dot model 

Our model has an exact correspondence with an ef- 
fective wide band two level Fano- Anderson model 2 ^ with 
Hamiltonian 



H 



e L d\d L 



■E 



£R.d R d R + T c d\d R + T^d R c 



Cka 



(13) 



k . a 



with a left and a right dot state, g?l and d R , and a cou- 
pling Tq between them. Each lead couples only to its 
adjacent dot. The transmission probability is 3 - 9 



r L r R \f c \ r < 



\(uj-e L + ir L /2)(u - e R + it R /2) - \T c f 



(14) 

Comparison to Eq. ([5|), with appropriate T R (cu) and 
A(uj), reveals the exact mapping T (uj) = T(lu) by the cor- 
respondence 



\fc\ 2 
e R 



r 2 

1 R 

2S R , 

£R- 



,o/(2tt), 



(15) 



III. NON-MARKOVIAN DYNAMICS 

A. Non-Markovian Master Equation (NMME) 

To treat non-Markovian effects in a transport master 
equation framework one writes the master equation as an 
integro-differential equation for the n-resolved reduced 
density matrix 

t 

PnW = E / W n-n'(t ~ ?) p n , (t')dt', (16) 



where n denotes the number of charges that have crossed 
the considered system. In Appendix [A] wc show how 
our system can be described with such an equation using 
the Born approximation. All W n (t) vanish except when 
n = or n = 1, where we find 



with 



Wo(t) = 
Wi(t) = 



-7l(«) 

lL(t) -l R {t) 

7fl(*) 




(17) 



2^ 



ri?(w)cos[(tj - ed)i] 



-^° e -fe*cos[( £i? -e d )t] 



7t(*) = 2 / — r iC os[(w-e d )t] 

= 2r L j(t). 



(18) 



We proceed with performing a Fourier summation and a 
Laplace transform via 



oo 

W(X,z) = ]TV™*| , 



'Wn(t)di 



W n (z) + Wi(z)e JX 

-7i(«) e lx 7fl,(2) 
Jl(z) -J R (z) 



(19) 



Since the cumulant generating function can be expressed 

solely in terms of the transmission probabilities ^ 40 ' 41 , not where 

only the first two cumulants 4 ^, but all current cumulants, 

of the two models coincide. 

When we plot current and Fano factor F = SJ (Ie) as 
a function of the detuning e<j — e R in Fig. [2] we find the 
typical structures of the corresponding quantities for the 
non-interacting double quantum do t 42 ' 43 : The current 
exhibits a maximum when the detuning vanishes. For 
the noise we find a minimum in resonance for broad tun- 
neling rates, when d R becomes sufficiently small a new 
maximum in the Fano factor at resonance grows out of 
the minimum. 



1r{z) 
1l(z) 



r 2 

1 Rfi 



z + S R 



it (z + S R ) 2 + (e R -e d ) 2) 



(20) 



We will use the abbreviation j(z) — Jl{z) + *f R (z). In 
order to avoid the tedious inverse Laplace transform we 
can use a recently developed elegant method to evalu- 
ate current and nois o 11 ' 44 . As it has only applied a few 
times until no w 11 ' 44 ' 45 ' 46 we show the explicit calculation 
for our model in appendix [B] and obtain the following 
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FIG. 2: The left plot shows the right tunneling rate Tr(oj) as a function of the energy u. We further show how the current I 
and Fano factor F behave as a function of the detuning ed — sr. The bandwidth Sr takes the three values shown on the right. 
For all curves we have r_L = r_R,o. Results from exact Green's functions. 



formulae: 



S 



Ie- 



7(0) 

■7 2 (0) 
7l(0) 

no) 



2^(0)] 
■2f«(0)]}. 



(21) 



We easily obtain the occupation of the dot by 
evaluating the density matrix in Laplace space, 
p(z) = [z - W(z)]- 1 p(t = 0). To obtain the time- 
resolved dynamics, one has to perform the inverse 
Laplace transform (Bromwich integral) by collecting all 
the corresponding residues. 

Markovian master equation (MME) — A Markovian 
master equation follows from (|16p by the integration 
C„ = | °° W n (t)dt = W n (z = 0) and leads to 



Co 



-r L o 



T R (e d ) 




(22) 



When we use this to evaluate noise we end up with the 
same result as in Eq. (j2Tj) . but without the derivatives 
7^(0) and 7^(0). For the equilibrium density matrix we 
find p 00 (t = 00) = T R (e d )/(T L + T R (e d )) , which is iden- 
tical with the NMME result. The advantage of the addi- 
tional Markov approximation is that the positivity of the 
density matrix will be conserved as one obtains a Lind- 
blad type master equation^. The disadvantage is that 
some information about the shape of the tunneling rates 
is lost. 



then interpolates through the solutions p T (t) = e £T '*po 
at t = t. The coarse grained Liouvillian can be derived^ 
by matching the second order expansion of the formal 
solution in the interaction picture x(t) = U (t)x{0)U^ (t) 

(where U(t) = Tcxp j-i J* H SB (t')dt'} with the time 

ordering operator T) with the second order expansion of 
p T {t) = e c t p T (0) at time t = r. For our specific model, 
following reference^ we obtain 



2tt' 



-Tr(w) 



(23) 

where sine x = The Coarse Graining method 

combines the two advantages of the Born and the 
Born-Markov approximation: For finite times it is 
sensitive to the shape of the tunneling rates, and 
at the same time it preserves positivity since the 
C T are of Lindblad form. Due to the identity 
linv^oo t ■ sine 2 ( UJ ~^ T — 2tt5(uj — e d ) we see that for 

large times the coarse graining method yields the same 
steady state density matrix as the two other approxi- 
mations. Also for current and Fano factor we have re- 
produced the results of the MME. In contrast to fixed 
graining-time derivations of master equationsi^, DCG 
dynamically adapts the coarse-graining time with the 
physical time, which in the long-time limit yields the 
Born Markov secular approximation. Both approaches 
yield completely positive maps but may lead to different 
stationary states. 



B. Dynamical Coarse Graining (DCG) 

A second approach to quantum transport is the re- 
cently developed dynamical coarse graining method^. 
The coarse graining method is also a second order weak 
coupling approximation although it can be extended 
to higher orders. Instead of solving a single master 
equation, it solves a continuous set p T (t) = C T p T (t) and 



IV. DISCUSSION 



A. Current 



To get a first idea of the difference between the exact 
solution and the NMME it is instructive to take a look 
at the stationary current /. The stationary current is 
not sensitive to non-Markovian effect&ii, since it only de- 
pends on the steady state occupation, and thus all three 
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FIG. 3: Current / as a function of the width Sr on 
a logarithmic scale in absence (upper part) and presence 
(lower part) of detuning, exact (red, solid) and with the 
NMME=MME=DCG (green, dashed) with the parameters 

Inset: For the marked points at Sr/Tr^o = 0.1 we plot as 
functions of ui the left tunneling rate Fl, the dot spectral 
function A(u), and the right tunneling rate Tr(uj), with a 
red density plot of the transmission probability T(ui) in the 
background. The dashed green line marks the level of the dot 
state. The master equations depend only on FR(td), while 
the exact solution is sensitive to the detailed shape of A(oj) 
and T(lu). 



presented approximations coincide. The two plots in Fig. 
[3] show the current as a function of the width of the right 
tunneling rate, Sr. Wc observe that for large Sr (Marko- 
vian limit) the master equations meet the exact solution 
very well, while for small Sr enormous deviations occur. 
In the absence of detuning, ed — sr : = 0, the master equa- 
tions overestimate the current, but they underestimate it 
for sufficiently large detuning. As an inset in Fig. [3] we 
choose the value Sr/Trq = 0.1 and show how the spec- 
tral function of the dot state A(w), the right tunneling 
rate Tr(u) and the transmission coefficient T(lo) behave 
there. 

Without detuning the spectral function exhibits a min- 
imum at energies where the tunneling rate is maxi- 
mal. Thus the true quantum mechanical transmission 
is smaller than the master equation result, which de- 
pends exclusively on the value of the tunneling rate at 
the dot level, Tr^)- This leads to the underestimation 
of the current by the master equations when detuning is 
present, because Tr^) is very small then. The exact 
result, by contrast, takes into account that the spectral 
function exhibits a small side peak near the maximum of 
Tr(uj) that enhances the probability for tunneling. 



B. Fano factor 

In Fig. |4] we show results for the Fano factor 
F = S/(Ie) (exact, with NMME and with MME/DCG). 
When the condition 



is fulfilled, all formalisms agree. Furthermore, for con- 
stant tunneling rates, the NMME yields the Markovian 
result F — r i +r n^ d ^ 

The exact Fano factor shows one minimum for small 
detuning and two minima for large detuning, which origi- 
nates from the double peak structure of the spectral func- 
tion shown in Fig. [3J However, there is no simple quan- 
titative connection between the locations of the extrema 
in the Fano factor and in the spectral function. This is 
similar to the Fano factor as function of the detuning 
in Fig. [5J where for small width Sr a second minimum 
appears, that is again not connected to the spectral func- 
tion's properties in a simple way. 

The NMME reproduces one of the two minima, but 
not both, and it produces supcr-Poissonian noise, where 
it should not appear. If the detuning e c i — sr becomes 
very small, the NMME even overestimates the minimum 
so strongly that it yields an unphysical negative Fano 
factor. 

The MME is not everywhere close to the exact result, 
but it is, on average, closer than the NMME result, and 
per construction yields physical results (Lindblad form). 
For the DCG method we have found that it yields the 
same result as the MME. 
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FIG. 5: Time dependent occupation probability rid(t), exact 
solution and solution in the three presented approximations 
Parameters are €d = Fh,o and td — en = Tl = 0.ir.R,o- The 
width 5r takes the values denoted in each plot. 
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FIG. 6: Frequency ojq as a function of the detuning uj — ujr 
for the parameters Fl ~ O.lTufi and = Ff^o- 



C. Time resolved occupation probabilities 

One key advantage of the DCG approach results from 
its ability to preserve positivity and at the same time to 
be more sensitive to the shape of the tunneling rates than 
the MME^£. We therefore evaluate the time dependent 
occupation probability of the localized level with four dif- 
ferent methods: exact, with the NMME, with the MME, 



and with DCG. In Fig. [5]wc find that at sufficiently wide 
bands, i.e., Sr > Tr.q, all three approximations meet the 
exact result very well. 

For smaller 5r one recognizes that in the stationary 
limit, the three approximations coincide with each other 
but not with the exact solution. The most prominent 
feature of the exact occupation rid(t) are oscillations as a 
function of time t. The NMME approximation captures 
the oscillations as non-Markovian features of the reser- 
voirs, as it should, but it strongly overestimates them. If 
the steady state is sufficiently close to zero, the NMME 
can lead to negative occupation probabilities, like in the 
third plot. 

In order to better understand the exact result, we ex- 
plicitly perform the two integrations in Eq. ([9]), which 
leads to 



Tld(t) = 

r 

I L 



4e 2 ) 



Mr 2 +4e 2 

(oji - e R + id R )(u>l 



-iTl 



-ir, 



(wi -u 2 ){u)* -w|)(o;i - 
(uj 2 - e R + i5 R ){ul - e R 



R) gi(i<J*— Wi)t 



Mr) 



(uji - W 2 )(Wi 
-2F L Im ( K * €R 



-uj 2 )(u 2 -wl) 

i$R){u2 - tR - 



iSr) 



R) e (iu> -V/2)t 



(0Ji - W 2 )(w* - <4) (w 2 - W*) 

(25) 

with the abbreviations r := T l + 25r and e := — eR. 
The frequencies lo\ and lo 2 are the poles of the spectral 
function A(u>) with negative imaginary part, and 01% and 
u>\ are its other two poles. The first part of n(t) in Eq. 
(|25[l corresponds to the steady state. The next two lines 
describe how the system performs exponential decay to- 
wards this steady state, because the exponents are real 
and negative. The last line of Eq. ([25)1 is responsible for 
the oscillations, because its exponent contains the imag- 
inary part iu>ot, which is given by 



Re< 



-4( £ R-iS R )(e d ~iT L /2) + A 



F 2 
1 R,0 

2tt 



1/2- 



.(26) 



If we switch off the coupling to the left side and the width 
on the right, i.e., — > and Sr —> 0, we recover the frc- 
I ~ 

quency: loq = y e 2 + 4-^ of coherent oscillations in an 
isolated two-level system (remember the mapping to the 
double dot model in Eq. (fTS"]) ). In Fig. [6] one recog- 
nizes that in the absence of the detuning — er, the 
frequency ujq completely vanishes with a nonanalyticity 
where Sr/Tr^q is of the order of 1. This nonanalyt- 
icity appears, when the radicand in Eq. (|26|) changes 
sign (at ed = sr it is purely real). When detuning is 
present, the frequency u>o does not completely vanish for 
S R — > co, but the oscillations of rid(t) are damped with 
the rate F/2 = S R + T^/2, such that in the Markovian 
limit, S R — > co, no oscillations survive. 
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V. CONCLUSIONS 

The main message of this paper is that one must be 
careful with non-Mar kovian master equations (NMMEs). 
We know that what we have here called "Markovian mas- 
ter equations" (MMEs) may lead to incorrect results 
when non-Markovian effects are strong. However, we 
have found that in the non-Markovian regime NMMEs 
can be worse than MMEs. On the one hand, NMMEs 
can give quantitative errors such as occur in Fig. ([4]). On 
the other, the positivity of the density matrix is not in 
general conserved, which can lead to unphysical results 
like a negative Fano factor and occupation probability, 
as in Fig. ((4]) and (|5|). We emphasize that the failure 
of the NMME does not result from the new evaluation 
techniques we have usedii, but instead from the way the 
Born approximation is performed in the derivation of Eq. 
(|T6| . We do not want to discourage using of NMMEs in 
general, but each time they are applied one should rea- 
son carefully why one expects them to yield better results 
than MMEs. 

We have tried the dynamical coarse graining (DCG) 
method as an alternative approach that is more sophis- 



ticated than the MME, but still in Lindblad form. For 
the time-independent quantities /, F and n(t — > oo) we 
have found no improvement in comparison to the MME. 
For the time-dependent occupation probability however, 
DCG does yield better results, at least for small times. 
We expect a similar improvement for the frequency- 
dependent Fano factor. 

For the model considered here, the condition <C 1 

' OR 

must be fulfilled to ensure a good quality of the three 
presented weak coupling approximations. It is certainly 
interesting to explore higher-order corrections to the pre- 
sented perturbative approximations. However, such ex- 
tensions should be treated with caution, since nice fea- 
tures such as positivity may be lost when higher orders 
arc included. 



Acknowledgments 

We would like to thank S. Gurvitz, C. Lopez-Moms, U. 
Klcinekathofer, and T. Novotny for helpful discussions. 
This work was supported by DFG grant BR 1528/5-1. 



* Electronic address: zedler@physik.tu-berlin.de 

1 Y. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 

2 Y. V. Nazarov, ed., Quantum Noise in Mesoscopic Physics 
(Kluwer Academic Publishers, Dordrecht, Boston, London, 
2003). 

3 Y. M. Blanter (2007), cond-mat/05 11478, to be published 
in Springer lecture notes, eds. Ch. Rothig, G. Schon and 
M. Vojta. 

4 G. KieBlich, E. Scholl, T. Brandes, F. Hohls, and R. J. 
Haug, Phys. Rev. Lett. 99, 206602 (2007). 

5 A. Braggio, C. Flindt, and T. Novotny, Physica E 40, 1745 
(2008). 

6 A. Braggio, J. Konig, and R. Fazio, Phys. Rev. Lett. 96, 
026805 (2006). 

7 I. Knezevic, Phys. Rev. B 77, 125301 (2008). 

8 R. Aguado and T. Brandes, Phys. Rev. Lett. 92, 206601 
(2004). 

9 T. Brandes, Phys. Rep. 408, 315 (2005). 

10 Y. Chen and G. Chen, Phys. Rev. B 77, 035312 (2008). 

11 C. Flindt, T. Novotny, A. Braggio, M. Sassetti, and A.-P. 
Jauho, Phys. Rev. Lett. 100, 150601 (2008). 

12 S. Datta, Electronic Transport in Mesoscopic Systems 
(Cambridge University Press, Cambridge, 1995). 

13 S. Datta, Quantum Transport: Atom to Transistor (Cam- 
bridge University Press, Cambridge, 2005). 

14 H. Haug and A.-P. Jauho, Quantum kinetics in transport 
and optics of semiconductors (Springer, Berlin, 2007). 

15 H.-P. Breuer and F. Petruccione, The theory of open quan- 
tum systems (Clarendon Press, Oxford, 2006). 

16 M. Nedjalkov, H. Kosina, S. Selberherr, C. Ringhofer, and 
D. K. Ferry, Phys. Rev. B 70, 115319 (2004). 

17 M. Nedjalkov, D. K. Ferry, C. Jacoboni, C. Ringhofer, 
I. Dimov, and V. Palankovski, Phys. Rev. B 74, 035311 
(2006). 



18 D. Querlioz, J. Saint-Martin, A. Bournel, and P. Dollfus, 
Phys. Rev. B 78, 165306 (2008). 

19 T. H. Stoof and Y. V. Nazarov, Phys. Rev. B 53, 1050 
(1996). 

20 S. A. Gurvitz and Y. S. Prager, Phys. Rev. B 53, 15932 
(1996). 

21 F. Kaiser, M. Strass, S. Kohler, and P. Hanggi, 
Chem. Phys. 322, 193 (2006). 

22 A. Thielmann, M. H. Hettler, J. Konig, and G. Schon, 
Phys. Rev. Lett. 95, 146806 (2005). 

23 X. Li, J. Luo, Y. Yang, P. Cui, and Y. Yan, Phys. Rev. B 
71, 205304 (2005). 

24 J. N. Pedersen and A. Wacker, Phys. Rev. B 72, 195330 
(2005). 

25 B. Elattari and S. A. Gurvitz, Phys. Rev. A 62, 032102 
(2000). 

26 W. Y. T. Matisse and W.-M. Zhang, Phys. Rev. B 78 
(2008). 

27 M. Thorwart, E. Paladino, and M. Grifoni, Chem. Phys. 
296, 333 (2004). 

28 P. Huang and H. Zheng, J. Phys. 20, 395233 (2008). 

29 U. Kleinekathofer, J. Chem. Phys. 121, 2505 (2004). 

30 P. W. Anderson, Phys. Rev. 124, 41 (1961). 

31 U. Fano, Phys. Rev. 124, 1866 (1961). 

32 G. D. Mahan, Many-Particle Physics (Springer Nether- 
lands, 2000). 

33 H. Bruus and K. Flensberg, Many-Body Quantum Theory 
in Condensed Matter Physics (Oxford University Press, 
2004). 

34 A. C. Hewson, The Kondo problem to heavy fermions 
(Cambridge University Press, Cambridge, 1993). 

35 A. D. Stone and P. A. Lee, Phys. Rev. Lett. 54, 1196 
(1985). 

36 M. Cini, Phys. Rev. B 22, 5887 (1980). 



8 



37 G. Stefanucci and C.-O. Almbladh, Phys. Rev. B 69, 
195318 (2004). 

38 G. Schaller, P. Zedler, and T. Brandes, Phys. Rev. A 79, 
032110 (2009). 

39 M. Y. Sumetskii and M. Fel'shtyn, JETP Lett. 53, 24 
(1991). 

40 L. S. Levitov and G. B. Lesovik, JETP Lett. 58, 230 
(1993). 

41 H. Lee, L. Levitov, and A. Y. Yakovets, Phys. Rev. B 51, 
4079 (1995). 

42 B. Elattari and S. A. Gurvitz, Phys. Lett. A 292, 289 
(2002). 

43 G. Kiefilich, P. Samuelsson, A. Wacker, and E. Scholl, 
Phys. Rev. B 73, 033312 (2006). 

44 C. Flindt, A. Braggio, and T. Novotny, AIP Conference 
Proceedings 922, 531 (2007). 

45 A. Braggio, C. Flindt, and T. Novotny, J. Stat. Mech. p. 
01048 (2009). 

46 D. Urban and J. Konig, Phys. Rev. B 79, 165319 (2009). 

47 G. Lindblad, Comm. Math. Phys. 48, 119 (1976). 

48 G. Schaller and T. Brandes, Phys. Rev. A 78, 022106 
(2008). 

49 D. A. Lidar, Z. Bihary, and K. B. Whaley, Chem. Phys. 
268, 35 (2001). 

APPENDIX A: DERIVATION OF THE KERNEL 

We decompose our Hamiltonian into the free 
part H Q = e d Sd + J2k,a e kacl a c ka and the coupling 

v = Y^k^kacft c ka + t*ka c \ a d)- This enables us to define 
the density matrix in the interaction picture, x(t). It 
obeys the Liouville-von-Neumann equation 



(Al) 



where V(i) = e lH aty e -iH t j g ^ e i n t erac tion picture ver- 
sion of the coupling. A tilde will mark the interaction pic- 
ture in all following text. By iterating Eq. (|A1|) twice, 
we obtain 



In the calculation that follows we find that the off- 
diagonal elements of p(t) decouple from the diagonal 
ones. With this knowledge we can choose the off- 
diagonals to be zero and neglect them in the density ma- 
trix, which we can therefore consider as a vector with two 
entries: 



(A5) 



( 



To evaluate the double commutator in equation (|A2|1 is 
lengthy but straight forward and yields 



pit) = J dt'W{t - t')p{t') 
o 



with 



W(t) 



-7i(*) mi*) 

lL{t) -lR{t) 



(A6) 



(A7) 



For the entries of W(i) we need to explicitly perform 
traces over the reservoirs: 

7i (t) = e fa **Tir{B2(0)Bi(t)iio} 

+e-^ t Tr{B 2 (t)B 1 (0)R }. (A8) 

With infinite bias, i.e., Rq = J^ fc <^ kL CkL this becomes 



(*) = ^|tfeL| 2 -2cos[(e fcL -e d )t] 



f ^T L (uj)cos[(u-e d )i\. (A9) 
Ztt 



In the same manner we get 

oc 



7fl(*) 



2 / -^T R (uj)cos[(uj-e d )t]. (A10) 
Ztt 



H(t) = -i[V(t),x(0)]- dt'[V(t),[V(t'),x(t')]}. (A2) 



At this point we perform the second-order weak coupling 
approximation by replacing the full density matrix 
by a tensor product of the reduced density matrix p(t) 
and the bath density matrix Rq, which we assume to be 
constant in time (Born approximation) . The partial trace 
over the first commutator in equation (|A2j) vanishes. To 
proceed we resolve the coupling into system operators Si 
and bath operators Bi such that 



where 
Si(t) 



V(t)=S 1 (t)B 1 (t) + B 2 (t)S 2 (t). 



tkaCka(t) : 

k,a 

d(t), B 2 (t) = ]TiL4aW- 

fc,a 



(A3) 



(A4) 



Physically, the upper right matrix clement of W(t) de- 
scribes a jump from the dot to the right lead. Thus 
we can say that it increases the number of passed elec- 
trons by one, while all other matrix elements leave it 
unchanged. Thus we distinguish 



W (i) 



-7l(«) 
idt) -Mt) J ' 



ma) = (° ?f> 



(AH) 
(A12) 



APPENDIX B: EVALUATION OF CURRENT 
AND NOISE 

We use a bra-ket like notation where the equilibrium 
state is represented by the ket 



|0»= limp(t) = J-( H°] 
" t^oo Hy ' 7(0) V 7i(0) 



(Bl) 



9 



that we obtain by setting p = in Eq. (fT6|) . We define 
the bra 



«0| = (1,1), 
and construct the projector 

V - V 2 - I0M0I - — ( ^ (0) ^ (0) 



(B2) 



(B3) 



and the projector Q = Q 2 = 1 — V . The resolvent of 
the Laplace-transformed kernel from Eq. (|19|) is 
K(A,x,z) := Q[W(x, z) - A ■ 1] _1 Q. We only need its 
value at zero, 1Z := 71(0, 0, 0), which is 



K 



1 



-7l(0) 7r(0) 



(B4) 



(7L(0)+7i?X0)) 2 V 7i(0) -7fl(0) 
We introduce the coefficients of the kernel's Taylor scries 



The expressions for current and the zero frequency 
noise 11 ' 44 are now obtained as 



I = e«6|W'|0»/z, 



((6|W'|o)}-2((6|WftW|o» A 2 

2ile \ (( 0|W'|0 )) - (( 0\WnW\0 )} 1 . (B6) 



W(x,z) = W + W'x + Wz 



2 ( W "X 2 + 2W'xz + Wz 2 ) + • • • (B5) 



Performing the derivatives and matrix multiplications, 
leads to the results of Eq. (T2~Tj) . 



